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Abstract. We study an influence of precise data on uncertainty of polarized parton distribution 
functions. This analysis includes the SLAC-E155 proton target data which are precise measure¬ 
ments. Polarized PDF uncertainties are estimated by using the Hessian matrix. We examine cor¬ 
relation effect between the antiquark and gluon uncertainties. It suggests that reducing the gluon 
uncertainty is needed to determine the polarized antiquark distribution clearly. 


INTRODUCTION 

Polarized parton distribution functions (polarized PDF’s) have so far been optimized 
from polarized deep inelastic scattering (polarized DIS) world data [|I], We could 
obtain only a slight piece of information about polarized antiquark and gluon distribu¬ 
tions. At this stage, the antiquark SU(3)^ flavor symmetry is assumed in most of the 
polarized PDF analyses. The SU(3)y^ symmetry breaking is already known as the Got¬ 
tfried sum rule violation in the unpolarized case. In principle, the polarized PDF analysis 
should take account of the symmetry breaking. However, we must not only determine 
a shape but also a sign of each polarized antiquark distribution. It needs more precise 
data to improve the current status. Semi-inclusive DIS experiments [@] are also expected 
to separate antiquark flavor distributions. However, the separated distributions may not 
be credible due to ambiguity of the fragmentation functions. Then, antiquark flavor dis¬ 
tributions cannot be decomposed clearly. The current knowledge of the polarized gluon 
distribution is still poor. The polarized gluon distribution is suggested as the positive 
distribution; however, there is large difference between various parameterization results. 

We would like to know ambiguity of polarized PDF’s quantitatively. PDF uncertainty 
plays an important role in illustrating the ambiguity. Furthermore, it is important to 
show the phenomenological uncertainty of predicted physical quantities (e.g., scattering 
cross-sections and spin asymmetries) with parameterized PDF’s and their uncertainties 
in our work. A purpose of this analysis is to clarify the current knowledge about the 
polarized PDF’s from the polarized DIS world data by using their PDF uncertainty. In 
this analysis, the polarized PDF’s are optimized including precise SLAC-E155 proton 
target data [^. Then, we examine an influence of the precise data on the polarized PDF 
uncertainty, which is estimated by the Hessian method. 


A fortran program of the AAC PDF library is available from |http://spin.riken.bnl.gov/aac; . 




PARAMETERIZATION OE THE POLARIZED PDF’S 


The polarized PDF’s are determined by using spin asymmetry of the polarized 
DIS experiments from the EMC, SMC, SLAC-E130, E142, E143, E154, E155, and 


HERMES: 


Ai(x,e2) 


2x[l+R{x,Q^)] 

F2ix,Q^) 


gi{x, 2 ^), 


( 1 ) 


where F 2 is the unpolarized structure function. The function R{x, Q^) = ojOj is the 
ratio of absorption cross sections for longitudinal and transverse virtual photons, and it 
is determined from experimental data in reasonably wide 2^ and x ranges in the SEAC 
experiments [Q]. The polarized structure function gj is expressed with polarized PDE’s: 
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where e- is the electric charge of quarks, and AC^, ACg are Wilson’s coefficient functions. 
The convolution 0 is defined by f{x) 0 g(x) = dy/y f{xly)g{y). The polarized PDE’s 
A/(= /^ — /^) are defined as helicity distributions in the nucleon. In the AAC analysis, 
the polarized PDE A/(x) is defined at initial 2^ by the weight function form: 

A/(x) = ( 1 + lx^)f{x) , (3) 

where f{x) is the unpolarized PDE, and A, a. A, and 7 are free parameters. Optimized 
PDE’s are four distributions; Auv{x), Adv{x), Aq{x), and Ag(x), and these are evolved 
from the initial 2^(=1 GeV^) to the same of experimental data by the DGEAP 
equation [0]. In particular, the gluon distribution Ag{x) contributes to the structure 
function with the non-zero coefficient function ACg in the NEO case. 

This analysis uses two constraint conditions. Eirst, the positivity condition is used 
to restrict large-x behavior of the polarized PDE’s. This condition corresponds to the 
probabilistic interpretation of the parton distributions in the EO: |A/(x)| < fix). It 
needs not to be satisfied strictly in the NEO analysis. However, the polarized antiquark 
and gluon distributions tend to badly break the positivity limit: |A/(x)| S> /(x). Such 
excessive behavior is due to the large experimental errors in the large-x region. Hence, 
this behavior should be limited by this condition. 

Next, the SU(3)^ flavor symmetry is assumed: Au{x) = Ad{x) = A 5 (x) = A 5 (x). Using 
this assumption, one can fix the first moments of the valence quarks with hyperon decay 
constants, then Auv = 0.926 and Adv — —0.341 are obtained. Note that the Bjorken sum 
rule is satisfied automatically by fixing first moments. Eurthermore, the spin content 
AE is obtained by AE^ ^3 = Awy + Adv -|- 6Aq. Since, the antiquark contribution is 

emphasized, then the spin content determination is susceptible to the antiquark behavior. 

In the analysis, we choose the modified minimal subtraction (MS) scheme, and the 
GRV parameterization for the unpolarized PDE’s at the NEO analysis [^. The total 
is minimized by the CERN subroutine MINUIT. 



UNCERTAINTY ESTIMATION 


Fortunately PDF uncertainty estimation method has been developed in the last several 
years (see a brief review [^). The polarized PDF uncertainty comes from several error 
sources, e.g., experimental errors, unpolarized PDF, and so on. However, it is 

difficult to incorporate these errors into uncertainty estimation simultaneously. In the 
present analysis, the polarized PDF uncertainty is estimated from experimental errors 
by using the Hessian matrix which is defined as a second order derivative matrix in 

the expanded function around its minimum point. The PDF uncertainty 5Af{x) 

can be obtained easily by the inverse matrix of the Hessian and linear error propagation: 
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where Ax^{= X^i‘^i) ~ Xmin) defined as the difference from the minimum x^- It deter¬ 
mines a confidence level of the PDF uncertainty, and it depends on the X^ distribution 
K{s) with N degrees of freedom. Here, N is the number of optimized parameters. In our 

estimation, the value of Ax^ is obtained by the following equation: K{s) ds = o , 

where 0 '(= 0.683) corresponds to 1 c error of a standard distribution in order to com¬ 
pare with general experimental errors. The statistical and systematic errors are added in 
quadrature, so that it could be overestimation. The proper estimation exists between the 
overestimated uncertainty and the uncertainty from only the statistical error. 


RESULTS AND DISCUSSIONS 

The best fitting result is X^i/d-o.f.) = 346.33(0.90). The first moments of new results 
and the AAC pervious results (AACOO, NLO set2) [|||] are shown in Table. [T|. A correla¬ 
tion coefficient between the first moment of the antiquark and gluon distributions is 
pqg — —0.836, and there is strong correlation between two distributions. The uncertain¬ 
ties of the new results become smaller than those of the previous results. The gluon first 
moment and spin content AE still have large uncertainty. The fixed first moments Auy and 
Ady do not have uncertainty, then the AE uncertainty is six times as large as the antiquark 
uncertainty. Thus, the spin content is subject to the uncertainty of the antiquark distri¬ 
bution. Figure |T] shows the uncertainty of the new antiquark distribution. The antiquark 
uncertainty becomes rather large in the region x < 0.01, however the experimental data 
scarcely exist. The polarized DIS spin asymmetries A^’'^(x) approaches rapidly to zero 
in the rang x < 0.004. It is insufficient to clarify small-x behavior of the antiquark dis¬ 
tribution. Therefore, the antiquark determination has extrapolating ambiguity in small-x 
region. It is needed tight constraint condition or other experiment. 

In addition. Figure |I] shows comparison between the PDF uncertainties of new results 
and the previous results. There are no significant improvements of the valence quark 
uncertainties. On the SU(3)y symmetry assumption, the fixing first moments strongly 
restricts the behavior of valence quark distributions. In contrast, the antiquark and gluon 
uncertainties are reduced in the range 0.01 < x < 0.5, where the E155 proton data exist. 




TABLE 1. First moments of the polarized antiquark, gluon, and 
spin content AL with their uncertainties at =1 GeV^. 



Aq 


AE 

New 

-0.062 ±0.023, 

0.499 ±1.268, 

0.213±0.138 

AACOO 

-0.057 ±0.038, 

0.532 ±1.949, 

0.241 ±0.228 


The precise polarized DIS data can reduce the antiquark uncertainty mainly. On the other 
hand, the gluon uncertainty changes in response to antiquark uncertainty reduction due 
to a strong correlation between two distributions. Since the gluon contribution to the 
structure function gj(x) is smaller than the quark and antiquark contributions, we can 
extract only a little information of the gluon distribution in spite of the NLO analysis. 
Actually, the gluon uncertainty is still large. It indicates the difficulty of determining the 
gluon distributions from the polarized DIS data. Therefore, the uncertainty reduction of 
the gluon distribution is due to the strong correlation rather than the NLO contribution. 

In order to examine the correlation effect on the parameterization, we re-analyzed 
the Ag(x) = 0 case in which the fixed gluon distribution does not have uncertainty. The 
polarized PDF uncertainties of the Ag(ji;) = 0 case are compared to those of the Ag(ji;) ^ 0 
case in Figure The gluon distribution slightly exists at high-2^ due to evolution 
of the singlet type DGLAP equation. The valence quark uncertainties scarcely change. 
Drastic improvement of the antiquark uncertainty is due to vanished the large gluon 
uncertainty, the obscure gluon distribution brings about the larger antiquark uncertainty 
by the complementary relation. 





FIGURE 1. Polarized PDF’s with their uncertainties at = 1 GeV^. Dashed curves are the uncertain¬ 
ties of previous results (AAC NLO-2) 




























FIGURE 2. Polarized antiquark and gluon distributions with their uncertainties at = 1 GeV^. The 
shaded portion shows the uncertainty of Ag(;c) = 0 results, and the dashed curves are the uncertainties of 
new results (Ag(x) ^ 0). 


SUMMARY 

By this analysis, the polarized PDF’s were optimized from the polarized DIS world 
data whieh ineluded the SLAC-E155 proton target data. The polarized PDF uneertain- 
ties were estimated by the Hessian method. The E155 preeise measurements seareely 
improve the valenee quark uneertainties, but they ean reduee the antiquark and gluon 
uneertainties. These, however, are still wrapped in large uneertainty. The SU(3)y^ sym¬ 
metry, whieh we are obliged to assume, restriets strongly the valenee quark behavior 
by fixing first moments, and the spin eontent determination depends on the antiquark 
behavior. Additionally, there is the strong eorrelation between the antiquark and gluon 
distributions. If the gluon distribution is elarified by RHIC-Spin at BNL, the uneertainty 
of the antiquark distribution ean be redueed to some extent. Similarly, the eomplemen- 
tary relation ean reduee the large uneertainty of the spin eontent whieh eomes from the 
extrapolating issue of the antiquark behavior. Then, we will be able to investigate the 
antiquark flavor dependenee in detail. 
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